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ABSTRACT 

We report on the measurement of the two-point correlation function and the pairwise 
peculiar velocity of galaxies in the IRAS PSCz survey. We compute these statistics first 
in redshift space, and then obtain the projected functions which have simple relations 
to the real-space correlation functions on the basis of the method developed earlier in 
analyzing the Las Campanas Redshift Survey (LCRS) by Jing, Mo, & Borner (1998). 
We find that the real space two-point correlation function can be fitted to a power law 
£(r) = (ro/r) 7 with 7 = 1.69 and ro = 3.70 /t Mpc. The pairwise peculiar velocity 
dispersion 012 (fp) is close to 400 km s" 1 at r p = 3h Mpc and decreases to about 
150kms _1 at r p ~ 0.2 /i _1 Mpc. These values are significantly lower than those obtained 
from the LCRS. 

In order to understand the implications of those measurements on the galaxy bias- 
ing, we construct mock samples for a low density spatially-flat cold dark matter model 
(r?o = 0.3, Ao = 0.7, T = 0.2, <7g = 1) using a set of high-resolution N-body simula- 
tions in a boxsize of 100/i _1 Mpc, 300/i _1 Mpc, and SOOft^Mpc. Applying a stronger 
cluster-underweight biasing (oc M -0,25 ) than for the LCRS (oc M -0,08 ), we are able 



- 2 - 



to reproduce these observational data, except for the strong decrease of the pairwise 
peculiar velocity at small separations. This is qualitatively ascribed to the different 
morphological mixture of galaxies in the two catalogues. Disk-dominated galaxy sam- 
ples drawn from the theoretically constructed GIF catalog yield results rather similar 
to our mock samples with the simple cluster-underweight biasing. We further apply the 
phenomenological biasing model in our N-body mock samples which takes account of 
the density-morphology relation of galaxies in clusters. The model does not reduce the 
velocity dispersions of galaxies to the level measured in the PSCz data either. Thus 
we conclude that the peculiar velocity dispersions of the PSCz galaxies require a bias- 
ing model which substantially reduces the peculiar velocity dispersion on small scales 
relative to their spatial clustering. 

Subject headings: galaxies: clustering - galaxies: distances and redshifts - large-scale 
structure of Universe - cosmology: theory - dark matter 

1. Introduction 

Large catalogs of galaxies with their redshift and angular positions are the major astronomical 
data sets which provide quantitative information on the distribution and formation of the universe. 
To properly decipher this information requires the use of statistical tools, and to assess their 
importance heavily relies on a quantitative comparison with theoretical models. 

In a recent series of papers (Jing, Mo, & Borner 1998; Jing & Borner 1998, 2001), various sta- 
tistical quantities from the Las Campanas Redshift Survey (LCRS) have been determined including 
the two-point correlation function (2PCF), £(r), the power spectrum, P(k), the pairwise peculiar 
velocity dispersion (PVD), ai2(r), and the three-point correlation function (3PCF), C( r i2, ^23, »"3i)- 
In addition, a detailed comparison between these observational results and the predictions of cur- 
rent cold dark matter (CDM) models has been carried out. Jing, Mo, & Borner (1998; JMB98 
hereafter) have constructed 60 mock samples for each theoretical model from a large set of high- 
resolution N-body simulations. The observational selection effects can then be taken into account 
in the analyses of the theoretical models exactly as for the real observational data. JMB98 have 
demonstrated at length that such a procedure is essential for the proper comparison between models 
and observations. 

JMB98 found that both the real-space 2PCF and the PVD can be measured reliably from the 
LCRS, and that the observed 2PCF for LCRS is significantly flatter than the mass 2PCF in CDM 
models on scales < l/t _1 Mpc. The observed PVD also turned out to be lower than that of the 
dark matter particles in these models. JMB98 introduced a cluster-underweight (CLW) biasing 
model to account for these discrepancies; CLW essentially assumes that the number of galaxies per 
unit dark matter mass in a massive halo of mass M decreases as oc M~ a . If a = 0.08, the 2PCF 
and the PVD of the LCRS are well reproduced in a spatially-flat CDM model with f^Vs ~ 0.4, 
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where is the density parameter of the model and as is the current rms linear density fluctuation 
within a sphere of radius 8/i~ 1 Mpc. In fact, Carlberg et al. (1996) argue that this anti-bias is at 
the level observed for rich clusters of galaxies. We note here that this CLW has been applied also 
in recent analytic modeling of the galaxy distribution (Seljak 2000; Peacock &; Smith 2000; Sheth 
et al. 2000). 

In the present paper, we report on the measurements of the 2PCF and PVD of galaxies 
in another large redshift survey, the PSCz catalog, which has become publicly available recently 
(Saunders et al. 2000). Since galaxies in the PSCz survey are selected from the IRAS point source 
catalog, they are supposed to be preferentially dominated by late types in contrast to the LCRS 
galaxies. Thus the difference of those clustering statistics between LRCS and PSCz should be 
interpreted as an indication for the morphology-dependent biasing of galaxies. 

The rest of the paper is organized as follows; we briefly describe the PSCz sample and the 
corresponding mock samples from simulations in §2. Then we present the results for the 2PCF of 
the PSCz in §3 which are compared with various model predictions in §4. Also shown there are the 
results for spiral galaxies predicted by the GIF simulation (Kauffmann et al. 1999a,b). Similarly 
we present in §5 and §6 the observational determination of PVD and its theoretical implications, 
respectively. Finally §7 contains a summary and further discussion. 

2. The PSCz galaxy sample and our mock catalogs 

The PSCz survey (Saunders et al. 2000) contains the redshifts for galaxies selected from the 
IRAS Point Source Catalog (PSC; Beichman et al. 1988). In total the PSCz catalogue contains 
15411 galaxies with their flux at 60/im exceeding 0.6Jy across 84% of the entire sky. The high sky 
coverage and the uniformity make this publicly available survey an ideal data set for the statistical 
studies of galaxy clustering. Both those parts of the sky excluded from the survey and the selection 
function are well documented in Saunders et al. (2000). In the analysis below, we use a volume- 
limited sample of galaxies with redshifts less than 6000 km s -1 according to the flux at 60/im and 
also all the galaxies with redshifts between 6000 km s -1 and 20000 km s _1 , resulting in a final semi- 
volume limited sample of 9425 galaxies. The volume-limited sample within 6000 km s -1 should 
reduce the influence of the local supercluster on the statistics on small scales, thus minimizing 
possible bias from the local supercluster. The redshift selection function is the same as that given 
by Saunders et al. (2000) for 6000 < z < 20000 km s _1 , is a constant for z < 6000 km s -1 and zero 
for z > 20000 km s" 1 . 

We construct a large set of mock samples, which have the same survey geometry, selection 
function and sampling rate as the real PSCz galaxies, from high-resolution cosmological N-body 
simulations. We emphasize that this is a very important aspect of our analysis in asserting the 
statistical significance of the results. For this purpose, we mainly consider a spatially-flat model 
with (f^o, Aq, r, 0"g)= (0.3, 0.7, 0.2, 1.0) simulated in a comoving boxsize of 300/i _1 Mpc, where 
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Ao is the cosmological constant and T = tloh is the shape parameter of the CDM power spectrum 
(Bardeen et al. 1986). We use three different realizations for each model employing 256 3 ~ 17 
million particles (Jing & Suto 1998; Jing 1998). This set of parameters is now considered as the 
standard low-density CDM model. JMB98 have found that this model reproduces the 2PCF and 
the PVD of the LCRS reasonably well. An even better fit to the LCRS data was obtained for a 
model with S7o = 0.2. The changes for this model can be roughly seen by a simple scaling of the 
term og^o' 6 which is proportional to the amplitude of the PVD. For reference, we also consider a 
Ao = Standard CDM model simulated in a 300 /i _1 Mpc box as well as additional LCDM models in 
lOO/i^Mpc and 800/t _1 Mpc to make sure that our conclusions are free from the numerical artifacts. 
The parameters for those simulation models are summarized in Table 1. 

The mock samples are generated in the same way as described in JMB98, using the selection 
function from the PSCz (Saunders et al. 2000). We construct 20 mock samples for each realization, 
and thus altogether 60 mock samples for the LCDM model in a 300ft, _1 Mpc box . For each mock 
sample we apply the CLW bias (JMB98). The CLW model, by construction, gives less weight to 
the denser clusters than the dark matter model, and thus is especially appropriate for the PSCz 
which is supposed to be dominated by spiral galaxies preferentially avoiding the high-density cluster 
regions. To be more specific, we randomly select galaxies from the dark matter particles so that the 
number of galaxies per unit dark matter mass, N/M, becomes proportional to M~ a within a halo 
of mass M. We apply this anti-biasing scheme for all dark halos with M > 7 x 10 11 h~ 1 M Q . While 
the model with a = 0.08 reproduces the 2PCF and the PVD for the LCRS data, we find (see below) 
that the PSCz data require a = 0.25, indicating stronger anti-biasing than the optically selected 
galaxies. We also apply this same bias model to the SCDM when we compare the predictions of 
this model to the observations. 

In addition, we also construct 20 mock catalogs from the GIF data base (Kauffmann et al. 
1999a, b) for comparison to our models and to the data. The GIF simulation adopts the same 
cosmological model, but with a slightly lower normalization of the linear power spectrum a% = 0.9. 
The simulation box size is 141/i -1 Mpc. The most important feature of the GIF data is that a 
semi-analytic galaxy formation model is implemented in contrast to our biasing scheme entirely on 
the basis of the gravitational clustering. We select "spiral galaxies" from the GIF mock samples 
according to their bulge-to-total luminosity ratio B/T. About 80% of those mock galaxies at the 
low end of B/T are retained for our analysis, taking into consideration the fact that the population 
of ellipticals is about 1/4 that of spirals in the local Universe. 

It might be more appropriate to select "galaxies" according to their infrared fluxes for the 
present purpose. Fortunately, however, the clustering properties of IRAS galaxies do not sensitively 
depend on the infrared luminosity (Szapudi et al. 2000), nor do those of the GIF spirals seem to 
depend sensitively on the ratio B/T. Thus our results below are fairly robust with respect to small 
changes in the population of galaxies. 
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3. The two-point correlation functions for the PSCz galaxies 

Quite a number of large galaxy catalogs, both two-dimensional (angular), and three-dimensional 
(redshift) catalogs, have been used to determine the 2PCF, and this statistic has now been estab- 
lished quite well. This statistic has produced constraints on theoretical models which are sensitive 
to the cosmological parameters, the initial power spectrum of the mass density fluctuations, and 
the biasing of specific galaxy types with respect to the dark matter, among others. 

We use the following estimator to measure the 2PCF £ z in redshift space; 

, , _ RR(r p ,7r) x DD(r p ,7T) 
UTp >* } ~ DR(r p ,Tr) 2 h Uj 

where DD is the count of galaxy-galaxy pairs in the projected separation r p and radial separation 
7T bins, RR and DR are similar counts of pairs formed by two random points and by one galaxy 
and one random point, respectively. 

We first display the redshift correlation function £(s) in Figure 1, which is the circular average 
of £z(?p,7r) for a radius s = ^jr 2 + ir 2 . The error bars have been estimated by the bootstrap 
resampling technique using the formula of Mo, Jing, k, Borner (1992). The resulting 2PCF is well 
fitted by £(s) = (5 /i -1 Mpc/s) L2 for s < 8/i -1 Mpc, but falls faster than the power-law on larger 
separations in good agreement with previous results (Seaborne et al. 1999). Here we focus our 
discussion on the projected 2PCF and the PVD (Davis & Peebles 1983), not only because these 
two statistics have not been determined for the PSCz galaxies before as of the writing of this paper 
but also because they provide very useful independent information about the spatial and velocity 
distribution of IRAS galaxies which can be used to understand and test galaxy formation models. 

We estimate the projected 2PCF w{r p ) from the integrated redshift 2PCF as 1 

poo 

w ( r p)= £,z{r p ,Tr)dir = yV 2 (rp,7Ti)A7ri, (2) 
Jo i 

where the summation runs from tt\ = 0.5/i -1 Mpc to -ztso = 49.5ft.~ 1 Mpc with AiTi = l/i _1 Mpc. 
JMB98 have shown that this method yields an unbiased estimate of w{r p ) which is also related to 
the real space 2PCF simply as 

w(r p ) = £ \(^rl+y 2 )dy. (3) 

Figure 2 plot the projected 2PCF for the PSCz with error bars estimated again with the bootstrap 
resampling technique. Our tests with the mock samples indicated that a standard scatter among 
different mock samples is actually comparable to these bootstrap errors. So we will use the bootstrap 



1 Note that our current definition of w(r p ) is a factor of two smaller than that of Davis & Peebles (1983), for 
example. 
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error for the observational estimates and the mock sample scatter for the comparison of the models 
with the observations. 

If the real space 2PCF £(r) is approximated as a power-law of the form: 

£(r) = (r /r)\ (4) 
then the projected 2PCF w(r p ) reduces to 

w(r P ) = v^ r(7 r / ( 2 y ~^ /2) (ro/r p yr p , (5) 

with T(x) being the Gamma function. As Figure 2 indicates, the observed w(r p ) of the PCSz 
galaxies for r p < 20/i~ 1 Mpc is fitted to the above form reasonably well with 

r = 3.70/i _1 Mpc, 7 = 1.69. (6) 

The slope is shallower and the amplitude is lower than those for the LCRS (JMB98). The error- 
bars are, however, substantially larger than for the LCRS, indicating that this catalog is somewhat 
noisier. The reason might be that the intrinsic clustering is weaker in this catalog and also that the 
volume surveyed by the catalog is smaller and in fact sparsely sampled. There is some indication for 
a deviation from the power-law between 0.2 h~ 1 Mpc and 0.8 h~ l Mpc, but it is not highly significant 
when considering the error-bars. The above best-fit parameters are in very good agreement with 
the result of the 1.2Jy sample, 7 = 1.66 and r = 3.76/t _1 Mpc (Fisher et al. 1994a). 



4. Comparison to the two-point correlation functions from the mock samples 

Figure 3 compares the projected 2PCF w(r p ) for the PSCz galaxies with those for the mock 
samples. The CLW bias model with a = 0.25 reproduces the PSCz data quite well within the la 
error (dotted lines); even the wiggly structure of the PSCz w(r p ) below 0.8 /i _1 Mpc is recovered. The 
value of a larger than for the optically selected galaxies (e.g., a = 0.08 for the LCRS) is consistent 
with the observational density-morphology relation; spiral galaxies preferentially selected by the 
IRAS PSCz do not concentrate in dense cluster regions. Moreover it is interesting that our CLW bias 
model can account for this behavior in a quantitative manner, purely by an appropriate adjustment 
of the parameter a. The SCDM model, when the CLW bias is applied, appears to underpredict 
clustering on scales of about 10/i~ 1 Mpc, as many previous studies have already found. 

To follow up this point a bit more, we have also constructed mock catalogs from the GIF 
simulation data (Kauffmann et al. 1999a, b). From the simulated catalogue, we select those galaxies 
with AVbg = Vb — V g > 1, where V& and V g denote the V-band magnitudes of the bulge and the 
whole galaxy. This encompasses about 80 percent of the galaxies in the GIF catalog. The resulting 
w(r p ) again fits the observations quite well (Fig.3) for r p £ l/i~ 1 Mpc. While the amplitude of 
w(r p ) for the GIF simulation is a bit larger than the PSCz data for r p £ 1 h~ 1 Mpc, it still lies close 
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to the +1(7 error line of our CLW mock samples. So these semi-analytic models of galaxy formation 
which incorporate physical processes like star formation and supernova explosions in some global 
way have a similar effect of reducing the number of galaxies per unit dark matter mass as our 
simple bias prescription. 

Figure 4 plots the 2PCFs in real space, £(r), of the GIF mock samples with different selection 
criteria, all of which follow a power-law oc r~ im quite nicely. While selecting galaxies on the basis 
of their AV& S reduces the 2PCF at r £ l/i _1 Mpc compared with those for all galaxies, the two 
different criteria, AVft 3 > 1 and AVf, g > 1.5, result in almost identical 2PCFs despite the fact that 
the latter case selects only 60% of the entire GIF galaxies. 

In Figure 5 we display several 2PCFs for dark matter from our simulations, the CLW mock 
samples with a = 0.08 and a = 0.25, and the AV& g > 1.0 sample from the GIF data set. Clearly, 
as has been noted many times already (see e.g. JMB98), the dark matter particles from the full 
simulation produce a £(r) which is too large and too steep for r < 1 h~ 1 M.pc. As for the LCRS 
galaxies, a scale-dependent bias is required for the LCDM model to be reconciled with the observed 
2PCF of the PSCz galaxies, which is consistent with the result obtained by Hamilton Sz Tegmark 
(2001) from an analysis of the real space power spectrum for the PSCz galaxies. The CLW bias 
model with a = 0.08 that fits the LCRS data well still predicts a much higher amplitude of £(r) 
than the PSCz data which are in good agreement with a = 0.25 model. The GIF spiral result lies 
between these two CLW bias models, and is closer to that of a = 0.25 at r ^ 1 ft^Mpc. We also 
note that the result of all GIF galaxies is quite close to the CLW model with a = 0.08. 

In summary, we may conclude that one can reproduce the 2PCFs both for the PSCz galaxies 
and for the optically selected LCRS galaxies in an LCDM model adopting a simple biasing pre- 
scription dependent on the type of galaxies in the sample, as long as the bias underweights galaxies 
in dense clusters. The semi-analytic GIF models may be a step on the way to provide a physical 
interpretation to the phenomenological CLW bias model. 



The Pairwise Velocity Dispersion (PVD) of galaxies is a well-defined statistical quantity which 
contains interesting information in principle on the cosmic matter distribution. The peculiar ve- 
locities of galaxies are determined by the action of the local gravitational fields, and thus they 
directly mirror the gravitational potentials caused by dark and luminous matter (Juszkiewicz et al. 
2000 and references therein). The PVD is measured by modeling the distortions in the observed 
redshift-space correlation function £ z (r p ,7r). The basic step in modeling is to write ^( r p; 7r ) as a 
convolution of the real-space 2PCF function £(r) and the distribution function fiyyi) of the relative 
peculiar velocity vyi of galaxy pairs along the line of sight (Davis & Peebles 1983): 



5. The pairwise peculiar velocity dispersion for the PSCz galaxies 
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While the real-space 2PCF £(r) can be directly estimated by inverting equation (3), we use the 
power-law fit (6) as in most previous work. For the velocity distribution function, we use the 
following exponential form which is supported from observations (Davis & Peebles 1983; Fisher et 
al. 1994b), theoretical models (Diaferio & Geller 1996; Sheth 1996; Seto & Yokoyama 1998), and 
direct simulations (Efstathiou et al. 1988; Magira, Jing, & Suto 1999): 



V2cr 12 \ Cl2 / 

where U\2 is the mean infall velocity and &\2 is the dispersion of the 1-D pairwise peculiar velocities 
along the line of sight. 

It is worth pointing out that the above modeling (eqs. [7] and [8]) is an approximation and 
the infall velocity vri is not known a priori, either. JMB98 demonstrated with the LCRS mock 
samples that the above procedure provides an accurate estimate of au (within 20% accuracy) if WpJ 
is known separately. In other words, vyi must be modeled carefully to allow a precise measurement 
of the dispersion a\2(r). In reality, however, U12 in the real Universe is not fairly certain at the 
present. One might think that on small scales UyJ is negligible, but this is true only for very small 
scales indeed. The function UyJ rises quite sharply around l/i _1 Mpc (Mo, Jing, & Borner 1997), 
reaching twice the Hubble velocity just beyond 1 /i _1 Mpc. Therefore the proper modeling of the 
infall velocity is crucial in an accurate measurement of the PVD au{r). 

Here we adopt the approach of JMB98 to determine the PVD for the PSCz catalogue. Figure 6 
displays the result of this modeling both for the case of a self-similar infall and for zero infall. The 
self-similar infall model has the same parameters as those used in JMB98, and is quite close to the 
mean pair velocity in CDM models with Signers = 0.4 for the relevant scales r < 5/i _1 Mpc. As 
emphasized in JMB98, only the use of some reasonable infall model (although the result is not very 
sensitive to the specific model used) gives a reliable reconstruction of the PVD, whereas the zero 
infall case does not even qualitatively reproduce the true value. We therefore consider the result 
of the self-similar infall as a reasonable estimate for the pairwise velocity dispersion of the PSCz 
galaxies. 

The procedure yields much lower PVD for the PSCz than that for the LCRS: oyi just reaches 
300 kms" 1 at r p = 1 h^Mpc, whereas <r 12 (l /i _1 Mpc) = 570±80km S - 1 for the LCRS (JMB98). 
This is again qualitatively consistent with the fact that spirals have smaller random motions than 
the galaxies that reside in big clusters. Our result at r p ~ 1 /i -1 Mpc is in good agreement with the 
value, 317 km s -1 , of Fisher et al. (1994b) for the 1.2Jy sample, although they did not show the 
scale-dependence of the PVD. 

Very remarkable in Figure 6 is the decrease of the PVD for small separations. The bootstrap 
errors in Figure 6 are quite substantial. Again we find that the errors are much larger than for the 
LCRS. Nevertheless it seems significant that the PVD decreases down to about 200kms _1 near 
r p = 0.2/i _1 Mpc. In contrast the LCRS data did not exhibit such a fall-off, but rather the PVD 
stayed almost constant near 400kms _1 . If we further take into account the observational error of 
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redshifts (typically ~ 120kms _1 ) and add it to the PVD measurement in quadrature, the real value 
becomes even smaller, ~ 150 km s -1 at r p = 0.2/i -1 Mpc. The galaxy pairs in the PSCz catalog are 
seemingly "very cold" ; they show very little random motion. 

6. Comparison with the pairwise velocity dispersion from the mock samples 

Figure 7 compares the PVD for the CLW and GIF mock samples with that for the PSCz 
galaxies. In all cases, the model predictions significantly exceed the estimate for the PSCz. Only 
around r p ps 3/i _1 Mpc the disagreement is not serious, especially if the large error bars are taken 
into account. In fact, we may speculate that a CDM model with S7o = 0.2 may even produce quite 
a good fit to the data at larger r p , since the amplitude of the PVD scales with crg^o' 6 - Reducing 
the values of o\2 accordingly, brings agreement on scales larger than 3fo _1 Mpc. There is, however, 
no way to reproduce the steep decrease towards small values (~ 150 km s -1 ) at r p = 0.2fo _1 Mpc 
of the PSCz data. The pairwise velocity dispersions are rather similar in the SCDM and in the 
LCDM models. 

This is really puzzling, because the LCRS galaxies do not show a signature like this. Even 
taking into account the large error bars at small scales, we have to concede a formal 3a devia- 
tion from the model results. This is consistent with previous findings of Ostriker & Suto (1990), 
Schlegel et al. (1994) and Strauss, Ostriker, & Cen (1998), who found that the galaxies in the local 
neighborhood are very quiet and that the dispersion of galaxies in mild density regions are small 
compared with the CDM predictions (also see Peebles 1995; Nolthenius, Klypin, & Primack 1997; 
Baker, Davis, & Lin 2000). It may also be that the measurements as well as the models become 
somewhat unreliable on such small scales. The finite sizes of galaxies may also reduce the PVD 
somewhat. As shown in Suto & Jing (1997) this can amount to up to 10 percent, if 0.1/i _1 Mpc 
is taken as a fiducial size, but the effect can be more pronounced for large and interacting spiral 
galaxies. 

Nonetheless our results present a new constraint on the well-defined PVD statistic that any 
successful galaxy formation theory should explain. At the moment, we may speculate that perhaps 
the spiral galaxies obey a very special velocity bias, or that the CDM models are way off in the 
description of small scale clustering as recently discussed in a different context (Spergel & Steinhardt 
2000; Yoshida et al. 2000). 

To separate clearly the effect of introducing the CLW bias and of constructing the mock 
catalogs, we present in Figures 8 and 9 the PVD from our full simulation for the CDM model and 
for the GIF simulation. The CDM models adopting the CLW bias with a = 0.08 and a = 0.25 
have a significantly lower PVD than the dark matter. We can also see that the shape of the PVD 
is quite similar among the dark matter, the CLW, and the GIF samples, while the amplitude of 
the dark matter particle PVD is typically higher by a few 100 km s" 1 . Again the GIF galaxies with 
AVbg > 1.0 or AVbg > 1.5 do not show any significant difference. The maximum in the PVD near 
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1 h Mpc present in the full simulation and the GIF sample is much reduced by the construction 
of the mock catalogs, and also for the a = 0.25 CLW models. 

So far we have used the 300 /i _1 Mpc simulations alone in the comparison. Since the boxsize is 
roughly the same with the PSCz survey itself once the outer boundary of the survey is concerned 
(in fact the weight of the most distant galaxies is small for our statistics and therefore the effective 
volume of the survey is much smaller than as the outer boundary indicates), our results may not 
faithfully reproduce the real cosmic variance even though we have three independent realizations. 
On the other hand, the small-scale peculiar velocity dispersions are supposed to be sensitive to 
the numerical resolution, and the 300/i _1 Mpc boxsize may be too large in this sense. To show 
that our 300/i _1 Mpc simulations are not statistically biased, we compute the projected 2PCF 
and the pairwise peculiar velocity dispersions for mock samples from 100/i~ 1 Mpc and 800/i _1 Mpc 
LCDM simulations also (again three realizations for each). The upper two panels of Figure 10 show 
the mean values of these two statistics. The agreement among the mock samples from different 
simulation boxes is satisfactory; the small difference between different simulations can be easily 
accounted for by the simulation parameters like the large softening length in the 800/t _1 Mpc. 
More important is the question if the error measurements from the 300ft, _1 Mpc simulations are 
underestimated. We show the ratio of the la error from the 800ft, _1 Mpc (100 /i _1 Mpc) simulations 
to that from the 300/i _1 Mpc simulations in the lower panels. The error ratios for both statistics 
and for both boxsizes are very close to 1 for the scales < 10/i _1 Mpc. This comparison clearly 
illustrates that the expected values and their sample-to-sample variations for w(r p ) and oyi are 
well represented in the three realizations in the 300/i _1 Mpc simulations and in the GIF mock 
samples. 



7. Discussion 

We have analyzed the data set of the IRAS PSCz galaxies (Saunders et al. 2000), and computed 
the two-point correlation functions and the pairwise peculiar velocity dispersion for these galaxies. 
A power-law fit to the real-space 2PCF £(r) = (ro/r) 7 gives an exponent 7 = 1.69, and a lower 
amplitude tq = 3.70/i~ 1 Mpc than the value obtained for the Las Campanas Redshift Survey (ro = 
5 /j^Mpc). We show that these results can be very well reproduced from mock samples constructed 
for a high resolution CDM simulation with f^o = 0.3 and Ao = 0.7 if we apply a biasing model which 
strongly underweights galaxies in the dense cluster regions. While this CLW bias was originally 
introduced by JMB98 to fit the LCRS galaxies, the bias needed for the PSCz galaxies is much 
stronger. More specifically, we find that the number of galaxies per unit dark matter in massive 
halos of mass M should be proportional to oc M~ a with a = 0.25 for the PCSz, as compared to 
a = 0.08 for the LCRS. The change of a from the LCRS to the PSCz is not unexpected, since 
IRAS galaxies tend to avoid high-density cluster regions. The biased mock catalogs fit the 2PCF 
of the PSCz galaxies extremely accurately. We also have selected samples from the GIF simulation 
and shown that appropriate data sets can be drawn from that numerical catalog which give similar 
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2PCF. 

The pairwise peculiar velocity dispersion measured from the PSCz has a much lower value, 
about 300 kms -1 at r p = l/i _1 Mpc, than the LCRS result at that separation of 570 ± SOkms" 1 . 
All the simulation models which are consistent with the 2PCF of the PSCz predict significantly 
larger values for PVD, although the CLW bias reduces the PVD to within the la limit of this value, 
at least near 3/i -1 Mpc (if f2o = 0.2). As discussed in the last section, the decrease of a\2{r p ) for 
r p % 1 h~ 1 Mpc for the PSCz data is significant, and cannot be reproduced by the models. To test 
the robustness of this important result, we applied another phenomenological biasing prescription 
to our simulations which mimics the observed density-morphology relation of galaxies in clusters. 
More specifically, we have selected particles within halos with the probability equal to the fraction 
of spiral galaxies in rich clusters (Dressier 1980). The fraction is set to be independent of the halo 
mass for simplicity. This biasing model does reduce both the 2PCF and PVD slightly, but the 
effect is much smaller than the difference between the CLW models and the PCSZ result. This 
strengthens our conclusion that the steep decline of PVD on small scales cannot be explained with 
the existing CDM model. 

In conclusion, we find that the infrared selected galaxies in the PSCz have a lower and less 
steep two-point correlation function than the optically selected sample of the LCRS. The PVD of 
the PSCz galaxies is also lower, and has a different r p dependence than the corresponding quantity 
measured from the LCRS. Spatially-flat CDM models with Qq ~ 0.3 can give a good fit to the 
data, except for the small-scale behavior of the PVD, if a strong cluster antibias is imposed. The 
infrared selected galaxies require a much stronger CLW bias than the optically selected ones. 
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referee for comments. Y.P.J, and G.B. are grateful for the hospitality extended toward them 
at Department of Physics, University of Tokyo, and for support during their stay at RESCEU. 
Y.P.J, is supported in part by the One-Hundred- Talent Program, by NKBRSF(G19990754) and by 
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Table 1. Simulation model parameters. 



Model 


n 


Ao 


rt 


0-8 


L hox [h x Mpc] 


r] gTav [h ^pc]* 


-^•particle 


"^particle [Mq] 


LCDM 


0.3 


0.7 


0.2 


1.0 


100 


39 


256 3 


5.0 X 10 9 


LCDM 


0.3 


0.7 


0.2 


1.0 


300 


117 


256 3 


1.3 x 10 11 


LCDM 


0.3 


0.7 


0.2 


1.0 


800 


400 


300 3 


1.5 x 10 12 


SCDM 


1.0 


0.0 


0.5 


0.6 


300 


117 


256 3 


4.5 x 10 11 


GIF 


0.3 


0.7 


0.21 


0.9 


141 


20 


256 3 


1.4 x 10 10 



t Shape parameter of the power spectrum, 
i Gravitational softening length. 
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Fig. 1. — The redshift two-point correlation function measured from the PSCz catalog (filled 
triangles). The straight line is a power-law (5 /i _1 Mpc/s) 1 - 2 . 
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Fig. 2. — The projected two-point correlation function measured from the PSCz catalog (filled 
triangles). Error bars are la deviations given by bootstrap resampling. The solid line is the best 
power-law fit. 
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Fig. 3. — The predictions of CDM models vs the observation for the projected two-point correlation 
function. Triangles show the observational result. The mean value and the la limits predicted by 
the cluster-weighted bias model are shown by the thick and thin lines respectively, and the mean 
values of the GIF simulation and the SCDM CLW model by the dot-dashed and dashed lines 
(without error bars). The SCDM curve is shifted vertically by a factor of 1/Vf to account for the 
necessary linear bias in this model. 
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Fig. 4. — The real space two-point correlation functions of the galaxies selected in the GIF simula- 
tion based on the magnitude difference AVb g between bulge and whole galaxy. The thin line is the 
fit f ( r ) = (3.7/r) L69 of the IRAS galaxies. 
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Fig. 5. — The real space two-point correlation functions of different bias tracers denned by the 
cluster-weighted model. For comparison, the result of the galaxies with AVb g > 1 in the GIF 
simulation are also plotted. The thin line is the fit = (3.7/r) L69 of the IRAS galaxies. 
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Fig. 6. — The pairwise velocity dispersion of the PSCz catalog determined for the self-similar infall 
(filled triangles) and for zero infall (open triangles). 
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Fig. 7. — The predictions of CDM models vs the observation for the pairwise velocity dispersion. 
Triangles show the observational result which has been corrected for the observational error of 
redshifts in quadrature. The mean value and the la limits predicted by the cluster-weighted bias 
model are shown by the thick and thin lines respectively, and the mean values of the GIF simulation 
and the SCDM CLW model by the dot-dashed and dashed lines (without error bars). 
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Fig. 8. — The pairwise velocity dispersion as determined from the 3-D velocity, of different bias 
tracers defined by the cluster-weighted model. For comparison, the result of the galaxies with 
AVbg > 1 in the GIF simulation are also plotted. 
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Fig. 9. — The pairwise velocity dispersion as determined from the 3-D velocity, of the galaxies 
selected in the GIF simulation based on the magnitude difference AV& 5 between bulge and whole 
galaxy. 
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Fig. 10. — Effect of boxsize in numerical simulations. 



